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ABSTRACT 

We report lifetime measurements of the 6 levels in the 3d*(^D)4d e^G term in Feii at an energy of 10.4 eV, and /-values for 14 
transitions from the investigated levels. The lifetimes were measured using time-resolved laser-induced fluorescence on ions in a 
laser-produced plasma. The high excitation energy, and the fact that the levels have the same parity as the the low-lying states directly 
populated in the plasma, necessitated the use of a two-photon excitation scheme. The probability for this process is greatly enhanced 
by the presence of the 3d®(^D)4p z®F levels at roughly half the energy difference. The /-values are obtained by combining the 
experimental lifetimes with branching fractions derived using relative intensities from a hollow cathode discharge lamp recorded with 
a Fourier transform spectrometer. The data is important for benchmarking atomic calculations of astrophysically important quantities 
and useful for spectroscopy of hot stars. 
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1. Introduction 


The present work reports on measurements of transition data in¬ 
volving high-excitation levels and lines of singly ionized iron, 
Fe II. The complex and line rich spectra of iron are observed in a 
wide variety of objects, such as stars, interstellar medium, solar 
coronae and quasars. In spectroscopy of stellar photospheres the 
atomic transitions most often appear as absorption lines. To use 
the spectral line for quantitative analysis of the population distri¬ 
bution and abundance determination in a stellar object, reliable 
/ values must be available. 

To meet the demands for accurate atomic data, ? initiated the 
FERRUM project with the goal of providing evaluated transi¬ 
tion data for astrophysical applications. Transitions from levels 
with different excitation potentials, and with a range of transi¬ 
tion probabilities have been measured and calculated for the iron 
group elements, see e.g. ???, and references therein. The present 
paper extends the data of highly excited levels, around 10 eV. 
The amount of data available from experimental studies is lim¬ 
ited, and the majority must be provided by calculations. The ex¬ 
perimental data is important to benchmark different calculations 
and to provide the user with reliable uncertainties and priorities 
between data sets. This is especially true for the complex ions 
where level mixing is important and hard to predict. 


A common technique to derive the oscillator strength is to 
combine measurements of the radiative lifetime of an upper 
level with relative intensities for all decay channels. The present 
paper reports on values using this approach, where the radia¬ 
tive lifetimes are determined with Time-Resolved Laser Induced 
Fluorescence (TR-FIF) from ions produced in a laser ablation 


plasma. These are combined with branching fractions obtained 
from intensity calibrated Fourier transform spectra using a hol¬ 
low cathode discharge as a light source. 

2. Radiative lifetime measurements 

Figure [2 presents a partial energy level diagram for Feii show¬ 
ing the large number of parent terms. The close similarity in the 
energy differences between 4s and 4p for all parent terms (the 
promotion energy), give rise to a high Fe ii line density in the 
region around 230-260 nm. In addition, the 4s-4p and 4p-5s/4d 
transition arrays have similar energy differences, and the latter 
transitions thus also appear in the same wavelength region. The 
target of the present study is the 3d®(^D)4d e®G term. This term 
has even parity which is the same as the ground conhguration 
3d®(^D)4s. This, and the large energy difference makes e®G im¬ 
possible to excite in a one-photon process. However, the use of 
an intense laser beam opens the possibility for two-photon exci¬ 
tations from the ground term, which has a high relative popula¬ 
tion in the ablation plasma. 

The TR-LIF set-up at Lund High Power Laser Facility is re¬ 
cently described in detail by ?, where a study similar to this, 
but for Cm, is reported; here we only emphasize the impor¬ 
tant aspects. The set-up contains two Nd:YAG lasers operating 
at 10 Hz. One of them (Continuum Surelite) is frequency dou¬ 
bled and focussed onto a rotating iron target placed inside a vac¬ 
uum chamber to generate the ablation plasma. The second laser 
is an injection seeded and Q-switched Continuum NY-82. The 
532 nm output from this laser was temporally compressed using 
stimulated Brillouin scattering in water before pumping a dye 
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laser (Continuum Nd-60) using DCIvf] To obtain the 237 nm 
radiation needed for the two-photon excitation, the output was 
frequency tripled in KH 2 PO 4 (KDP) and BaB 204 (BBO) crys¬ 
tals and Raman shifted in a high pressure hydrogen cell. The 
laser-induced fluorescence was detected by a 1/8 m monochro¬ 
mator with a 280 pm wide entrance slit oriented parallel to the 
excitation laser beam and perpendicular to the ablation laser. All 
measurements were performed in the second spectral order giv¬ 
ing an instrumental line width of 0.5 nm. The fluorescence signal 
was recorded with a micro-channel plate photomultiplier (PM) 
tube (Hamamatsu R3809U) with a rise time of 0.2 ns, and digi¬ 
tized by a Tektronix oscilloscope with 2.5 GHz analogue band¬ 
width. A second channel on the oscilloscope sampled the exci¬ 
tation laser pulse shape from a fast diode fed by scattered light 
from an attenuating filter. Each decay curve, as well as the tem¬ 
poral shape of the excitation pulse, was averaged over 1000 laser 
shots. The lifetimes were extracted using the program DECEIT 
(?). The program derives the lifetime from the measured laser 
pulse and the recorded decay. In the present experiment, the flu¬ 
orescence data was fitted with a single exponential convoluted 
with the square of the measured laser pulse (since the upper state 
is excited through a two-photon absorption) and a constant back¬ 
ground. The laser pulse had a temporal width (EWHM) of 2 ns. 
A typical example is shown in Figure!^ 

The final lifetimes given in Tablef^re the averages of around 
20 measurements, performed during different days. The quoted 
uncertainties, between 10 and 17%, include the statistical con¬ 
tributions from the fitting process as well as variations in the 
results between the different measurements. The high probabil¬ 
ity of spectral coincidences between the two-photon pumping 
channel and a single photon excitation of the 4p levels had to be 
thoroughly considered. Another consequence of the similarity 
between the 4s-4p and 4p-4d wavelengths is that the secondary 
cascades (4s-4p) coincide with the measured primary decay (4p- 
4d). The chosen pump and detection channels are given in Ta- 
bl^ and the most difficult cases will be discussed in detail be¬ 
low. All measurements are checked for blending lines, both in 
the excitation and decay channels, and three of the levels are dis¬ 
cussed in more detail below. 

The high J value of the e®Gi 3/2 level, and the AJ - 0, +1 re¬ 
striction for El transitions, effectively limits the decay channels 
to only the (^D)4p level with 7=11/2, z®Fii/ 2 . Transitions to 4p 
levels belonging to other parent terms are several orders of mag¬ 
nitudes weaker. The pumping transition is located at 238.0 nm, 
and the primary and secondary decays fall at 237.64 nm and 
238.28 nm, respectively. In this case the detection had to be 
shifted towards shorter wavelengths to avoid the influence from 
scattered laser light and secondary decay. However, shifting the 
detection wavelength too much to the blue brings in the fluores¬ 
cence at 234.8 nm from 4p z"^D 7/2 with a lifetime of 3.02 ns (?), 
which is excited by a single 238.0 nm photon. This effect was 
clearly observed in the decay curve as a weak multi-exponential 
shape which resulted in a longer lifetime in our single exponen¬ 
tial analysis compared to detection at longer wavelengths. By 
stepping the detection wavelength from this decay towards the 
primary decay at 237.64 nm, the contribution from z^D 7/2 could 
be monitored, both in the derived lifetime and the fitted curve. 
From this behaviour and the known instrumental line width of 
0.5 nm we estimated the maximum contribution from other de¬ 
cays to be 10%. This is reflected in the larger total uncertainty 
quoted for this level, 14%, compared to 11% for the other levels. 


* 4-(Dicyanomethylene)-2-methyl-6-(4-dimethylaminostyryl)-4H- 
pyran 


The level e®G 11/2 was pumped at 237.2 nm from the ground 
state a®D 9 / 2 . The primary decay is at 237.0 nm whereas the sec¬ 
ondary decay, a®D-z®F, falls at 238.2 nm. The monochromator 
was tuned to the blue side of the decay channel to minimize 
scattered laser light in the detection. In addition, the contribu¬ 
tion from secondary decays was avoided. 

The level e®G 7/2 can be populated by two-photon excitation 
from the lower levels a-^^ 3 / 2 , 5 / 2 j/ 2 - The excitation from 7=3/2 
was the best channel. The transitions from 7=7/2 is close to the 
excitation route to e®Gi 3 / 2 , but still possible. The transition from 
7=5/2 is close to the laser frequency and could not be used. 



Fig. 1. A partial level diagram of Fe ii showing the parent term struc¬ 
ture, and the similar energies for all 4s-4p and 4p-5s/4d transitions. The 
arrow marks the transitions investigated in this study. 


3. Branching fractions 

The BF for a transition from an upper level i to a lower level k 
can be expressed as: 


Tj Tj hj 

7=1 7=1 

Here N is the total number of transitions from level i and j 
is summed over all lower levels. The first equality is the defini¬ 
tion of BFik, i.e. the fraction of the atoms in level i that decay 
through the channel ik . The second equality assumes the case of 
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Table 1. Excitation and detection scheme for the measured levels, and results from this work and comparisons with previous calculations. 


Level 

E ! cm ^ 

Pump channel 

Level A / nm 

Detection channel 

A / nm 

Lifetime / ns 
This work 

K13" 

RU98* 

e'’Gi3/2 

84035 

a'’D9/2 

237.99 

237 

2.1+0.3 

1.48 

1.57 

e^Gii/2 

84296 

a®D9/2 

237.26 

236-237 

2 .0+ 0.2 

1.46 

1.56 

e®G9/2 

84527 

a®D7/2 

237.69 

235-238 

1.9+ 0.2 

1.45 

1.55 

e®G7/2 

84710 

a®D3/2,7/2 

237.10; 238.45 

235-263 

1 .8+0.2 

1.55 

1.55 

e^G5/2 

84844 

a®D5/2 

237.52 

235-266 

1 .8+0.2 

1.44 

1.54 

e®G3/2 

84938 

a®Di/2 

238.21 

235-236 

1.7+0.3 

1.43 

1.53 


a*} h*}*} 




Wavenumber / cm' 


Fig. 2. The first 20 ns of the decay of the 3d®4d e®Gii /2 level in Feii fol¬ 
lowing two-photon excitation from the ground state (3d^4s a®D 9 / 2 ). The 
evaluated lifetime is 2.0±0.2 ns. Background subtracted data points (-I-) 
are plotted together with a fitted single exponential decay convoluted 
by the square of the measured laser pulse (solid line). The dashed curve 
shows the squared laser pulse. 


Fig. 3. The 4p z®Fii /2 - 4d e®Gn /2 line at 42328.78 cm^'blended by two 
other Fell lines at 0.185 and 0.352 cm^'higher wavenumber, respec¬ 
tively. These separations have been kept fixed in the fitting of the three 
Gaussian shaped lines. 


1 


an optically thin plasma, where the observed photon intensity of 
a line, 7,^ measured as photons per second, is proportional to A,y 
(?) . When all intensities are measured on a common calibrated 
scale, BFii^ can be derived from the measured intensities. 

The lifetime of the upper level, t,-, puts the relative intensities 
on an absolute scale. With 

1 

- 

2 Au 

7=1 

and Eq. Q we derive A from 

A., = ^ ( 2 ) 

Ti 

This requires all lines from the upper level to be measured. In 
practice it is usually impossible to include all transitions, either 
because they occur outside the detector range or they are too 
weak to be observed. Thus, we rewrite Eq. Q as: 


BFit = 


lik 


N 

X lij 
7=1 


1 


+ 2 
— 


+ 


N 


X 

j=n+\ 


k 

lit 


( 3 ) 




N 


X Aij 

_/=w+l 


Here n is the number of observed lines and the sum over the 
intensity ratios for the unobserved lines is approximated using 
the theoretical transition probabilities from ?. 

The iron plasma was produced in a hollow cathode discharge 
(HCD), operated with neon as carrier gas and using currents be¬ 
tween 0.6 and 0.8 A and pressures between 1 and 2 torr. The 
light emitted was recorded with a Eourier transform spectrometer 
(Chelsea Instmments ET500) in the range 27640 - 52290 cm ' 
and detected by a solar blind PM tube (Hamamatsu R166). In 
total, five spectra were measured. While a higher current in the 
HCD produced a more intense spectrum it also led to slightly 
wider lines due to the temperature dependent Doppler effect. 
Line widths between 0.18 and 0.25 cm^^were observed with 
most of the spectra having 0.21 cm *. To calibrate the intensity 
scale, three spectra with lower resolution were recorded using a 
deuterium lamp with known relative spectral radiance, measured 
at the Physicalisch-Technische Bundesanstalt, Berlin, Germany. 

Despite the high spectral resolution of the spectrometer (re¬ 
solving power about 2 ■ 10^) the very line rich spectrum of Ee^ 
in this region produced a number of blends that must be taken 
into account to derive accurate experimental intensities. The first 
case is the z®Eii /2 - e®Gii /2 channel at 42328.78 cm“\ where 
there are two additional Ee ii lines at 0.185 and 0.352 cm 'higher 
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Wavelength (A) Wavelength (A) Wavelength (A) 

Fig. 5. Partial absorption spectrum of the chemically peculiar star chi Lupi, with the investigated Fe ii lines marked with tickmarks. The spectrum 
is observed by the Flubble Space Telescope and its Goddard Fligh Resolution Spectrograph, and is obtained from the MAST data archive at the 
Space Telescope Science Institute. See ? for a complete identification list and a comparison with a synthetic spectrum. 



Fig. 4. Experimental lifetimes are plotted with error bars (+), whereas 
the theoretical lifetimes (?) are marked with crosses (x). 


wavenumbers, respectively (??). This blending was handled by 
fitting a sum of three Gaussian shaped line components, con¬ 
strained to have the same FWHM and their separations fixed to 
the accurately known values, using the computer code GFit (?). 
One of these fits is shown in Figure]^ Consistent line widths and 
relative intensities were obtained in all five measurements. 

A more severe blending is found in the z^Fs /2 - e®G 7/2 
channel at 42375.86 cm 'where another Feii line is only 
0.11 cm 'higher. Furthermore, contrary to the previous case 
the blending line is also significantly more intense. The exper¬ 
imental spectra were analysed using the same technique of con¬ 
strained fitting as before. Also, in this case, we obtained the same 
FWFIM as for resolved lines but the relative intensities of the two 
components varied somewhat between the five measurements. In 
addition to the three lines included in our work, ? predicts an A 
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value larger than 10^ s“' for four additional transitions (41472, 
40478, 40263 and 39925 cm“') from the e®G 7 / 2 . Unfortunately 
all these lines are too severely blended in our spectra to be in¬ 
cluded in our analysis. This causes the so-called residual in Table 
2, i.e. the sum of all theoretical BFs from ? for the unobserved 
lines from this level, to be as high as 34%. Since the experimen¬ 
tal BFs, according to Eq (3), depends on this residual, our values 
for the e®G 7/2 level must be considered less reliable than for the 
other levels which is reflected in the higher uncertainties. 

The last example of blending is the close coincidence be¬ 
tween z®F 7/2 - e®G 5/2 at 42607.80 cm^'and a Feii transition at 
42607.79 cm '. In addition, the lines can be expected to have 
similar intensities. Thus the 7/2 - 5/2 combination are omitted 
from the BF analysis. 

We present the final average BFs, and transition probabilities 
obtained in Table where we include the theoretical results by 
? for comparison. The error estimate in the BFs, includes the un¬ 
certainty in the intensity calibration, in the area determination of 
the peaks, and a 50 % uncertainty in the contribution from the 
unobserved branches obtained from ?, as described in Eq (|^. 
Since all but one of the observed decay branches are within 735 
cm 'of each other and, in addition, occur in a wavenumber re¬ 
gion where the deuterium lamp calibration spectrum is quite flat, 
the uncertainty contribution from the calibration is small, about 
3 %. However, in the present case the total uncertainty of the 
transition probabilities for the strong lines is dominated by the 
contribution from the lifetime measurements. Eor the weak lines, 
the uncertainty in the BF is the dominating source. A more gen¬ 
eral discussion of uncertainties in experimental transition proba¬ 
bilities can be found in ?. 


4. Discussion 

Eigure|^shows a Hubble Space Telescope spectrum of the binary 
system chi Lupi (HD141556), which has spectral types B9.5p 
HgMn and A2 Vm with effective temperatures of 7’eff=10 650 K 
and 9200 K, respectively (?). The investigated transitions 4p- 
4d e®G are observed as prominent absorption features around 
235 nm. 

In Table we compare our experimental lifetimes with the 
available theoretical data from ? and ?. ? applied a superposition 
of configurations method using a modified version of the Cowan 
codes (?), whereas ? applied the orthogonal operator formalism. 
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Both calculations make use of experimental level energies to im¬ 
prove the results. As seen in Table [T] the two theoretical results 
are within 7%, but they are about 25 % shorter and do not over¬ 
lap with the experimental values, even taking the rather large 
uncertainties into account. A graphical comparison is shown in 
Figure]^ where a weak tendency of lower experimental lifetimes 
for lower J values is seen. This trend is within the error bars. 

The discrepancy between measured and calculated BFs is 
relatively small for most lines. An interesting exception are the 
transitions from e®G 7 / 2 , where there is a large difference between 
the two theoretical studies. Although both give similar lifetimes 
they differ significantly in the transition probabilities for indi¬ 
vidual lines. For two of the transitions our values support those 
of ?, whereas the z"^P 5/2 - 6 ^ 67/2 line at 37743 cm“', which is 
predicted to be quite intense is barely observable in our spec¬ 
tra. In addition, the large residual for this level, as discussed 
above, leads to a large uncertainty in the experimental A val¬ 
ues. The fact that the values by ? differ significantly from his 
earlier 2010 calculations indicate theoretical difficulties for this 
level. However, for all levels the resulting experimental A val¬ 
ues are a combination of lifetimes and BFs, and the discrepancy 
between the experimental and theoretical results is in most cases 
a consequence of the lifetime differences. 

The general agreement, except for the mentioned e®G 7 / 2 , of 
the experimental / values with the calculations by ? and ? pro¬ 
vides us with confidence to use either set for abundance deter¬ 
minations and also for high-excitation lines represented by the 
transitions and levels investigated in this work. We have a pref¬ 
erence for the calculations by ? when compared with the ex¬ 
perimental values, however, the discrepancy between experiment 
and theory has to be considered when estimating the uncertainty 
in the abundance determination or evaluation of the stellar at¬ 
mosphere models. We recommend that the experimental data is 
used where available, and complemented with theoretical calcu¬ 
lations for the other lines. 
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^ The calculations from 2010 (referenced as KIO in Kurucz notation) 
are not available online any longer, but replaced by the new data ?. 
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Table 2. The resulting transition probabilities and oscillator strengths for the lines studied. 
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® Uncertainty in the experimental A values 

f Residual is the sum of all theoretical BFs for the unobserved lines from this level. 






